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We discuss the possibility to generate, manipulate, and probe single spins in 
single-level quantum dots coupled to ferromagnetic leads. The spin-polarized 
currents flowing between dot and leads lead to a non-equilibrium spin accu- 
mulation, i.e., a finite polarization of the dot spin. Both the magnitude and 
the direction of the dot's spin polarization depends on the magnetic properties 
of leads and their coupling to the dot. They can be, furthermore, manipulated 
by either an externally applied magnetic field or an intrinsically present ex- 
change field that arises due to the tunnel coupling of the strongly-interacting 
quantum-dot states to spin-polarized leads. The exchange field can be tuned 
by both the gate and bias voltage, which, therefore, provide convenient han- 
dles to manipulate the quantum-dot spin. Since the transmission through the 
quantum-dot spin valve sensitively depends on the state of the quantum-dot 
spin, all the dynamics of the latter is reflected in the transport properties of 
the device. 



1 Introduction 

The study of single spins in quantum-dot spin valves resides in the inter- 
section of the two highly-interesting and extensively-pursued research fields 
of spintronics on the one hand and transport through nanostructures on the 
other hand side. Quantum dots consist of a small confined island with a low 
capacity such that a macroscopic gate or bias voltage is needed to add a single 
electron, leading to Coulomb-blockade phenomena 001 Ej . The notion that 
not only the charge but, simultaneously, also the spin degree of freedom of the 
electrons can be made use of, for example by using ferromagnetic leads, de- 
fines the field of spintronics 0J|S]. A quantum-dot spin valve, i.e., a quantum 
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dot coupled to ferromagnetic leads, exploits both the spin polarization of the 
electrons and the sensitivity of the charge to Coulomb interaction. Therefore 
electronic transport is governed by the behavior of a single spin. To discuss 
the possibility to generate, manipulate, and probe single spins via electronic 
transport through quantum-dot spin valves is the goal of this chapter. 

The capacity C of a metallic or semiconductor island decreases when 
shrinking its size. For small quantum dots, the energy scale to add a sin- 
gle electron on the island, the charging energy U = e 2 /2C, exceeds the energy 
scales set by temperature k-&T or bias voltage eV, and Coulomb-blockade phe- 
nomena arise, as first observed by Fulton and Dolan [Hj. If, in addition, the 
island size becomes comparable to the Fermi wavelength then the level spec- 
trum on the island will be discrete. For sufficiently large energy-level spacings, 
only a single level may participate in transport. Such as system can then be 
described by the Anderson-impurity model, which is introduced below. 

Famous examples of spintronics devices are the spin valves based on cither 
the giant magnetoresistance effect [7] in magnetic multilayers or the tunnel 
magnctorcsistance jH] in magnetic tunnel junctions. These effects arise, when 
two ferromagnetic leads are in contact via a conducting layer or a tunnel bar- 
rier, respectively. The transport characteristics of the device then depend on 
the relative orientations of the lead magnetizations. If, in a magnetic tun- 
nel junction, the lead magnetizations enclose the angle 4>, the conductance 
through the tunnel junction is proportional to cos^> |?flllO|. i.e., it is maximal 
for parallel and minimal for antiparallel alignment of the leads' magnetiza- 
tions. This angular dependence simply reflects the overlap of the spinor part 
of the majority-spin wave functions in the source and drain electrode, which 
is given by the externally controlled leads' magnetizations. 

This picture changes once spin accumulation can occur. Let us consider 
transport through a ferromagnet - nonmagnet - fcrromagnet sandwich struc- 
ture with the thickness of the normal layer being smaller than the spin diffu- 
sion length. A finite bias voltage applied between the two ferromagnets with 
nonparallel magnetization directions leads to a local imbalance of spin-up and 
spin-down electrons in the nonmagnetic layer. This non-equilibrium polariza- 
tion of the electrons in the nonmagnetic region, known as spin accumulation, 
mediates the information of the relative orientation of the leads' magnetiza- 
tion through the middle part, such that the transmission through the device 
is reduced for increasing angle <fr between directions of the leads' magnetic 
moments. 

An extreme limit of the above scheme is realized in a quantum-dot spin 
valve. It consists of a quantum dot that is tunnel coupled to ferromagnetic, see 
Fig. H In this case, the information about the relative leads' magnetization 
directions is mediated by a single quantum-dot electron that, as a conse- 
quence of a finite bias voltage is partially spin polarized, described by a finite 
quantum-statistical average S = (%/2)((t) of the dot spin. It is the orienta- 
tion of the dot spin relative to the leads together with the degree of the dot 
spin polarization that determines the transport, rather than just the relative 
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orientation of the leads' magnetization directions only. Any manipulation of 
the dot spin polarization will change the transmission through the device. 




Fig. 1. A quantum dot contacted by ferromagnetic leads with non-collinear mag- 
netizations. The lead magnetization directions enclose the angle 4>. By forcing a 
current through the system, a non-equilibrium spin S accumulates on the otherwise 
non-magnetic dot. 



A quantum-dot spin valve is, thus, a convenient tool to generate, manip- 
ulate, and detect spin polarization of single quantum-dot electrons. Both the 
generation and the detection of spin polarization on the quantum dot occur via 
electrical transport as a consequence of spin-polarized charge currents from 
and to the leads. One of the intriguing features of a quantum-dot spin-valve 
device is the possibility to further manipulate the dot spin. This can be done 
directly via an externally applied magnetic field But also the gate and 
transport voltages influence the dot spin [1211131 IHlll5| . To understand this, 
it is important to notice that the strong Coulomb interaction on the quan- 
tum dot yields many-body correlations. As for the spin degree of freedom, 
the quantum-dot electrons are subject to an exchange field that arises as a 
many-body effect due to the tunnel coupling to spin-polarized leads. This ex- 
change field sensitively depends on the system parameters such as the gate 
and bias voltage. The latter, therefore, provide suitable handles to manipulate 
the quantum-dot spin. 

Combining ferromagnetic (typically metallic) leads with quantum dots, 
which are usually semiconductor structures is experimentally challenging. Re- 
cent experimental approaches to such a quantum-dot spin valve involve metal- 
lic islands El > granular systems > carbon nanotubes p5] as well as 
single molecules [2U] or self-assembled quantum dots [21 HI] coupled to ferro- 
magnetic leads. Another possible realization would rely on contacting a surface 
impurity acting as quantum dot with a spin-polarized STM tip |23j . 

Successful demonstration of tunnel magnetoresistance through a strongly 
interacting system has been reported by Sahoo et al. [2I_ in single-wall car- 
bon nanotubes contacted by PdNi leads, and by Zhang et al. in Al grains 
sandwiched inside a tunnel junction between to Co leads [2*5| . 

The article is organized as follows: In Sec. we define the Hamiltonian of 
the quantum dot coupled to ferromagnetic leads. In Sec.|31we address the dy- 
namics of the dot spin and charge. Starting from a rigid calculation of the spin 
and charge current through a tunnel junction, we construct the master/Bloch 
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equation for the charge/spin degrees of freedom from the charge/spin con- 
tinuity equation. From the Bloch equation, we discuss, how to prepare, and 
modify the dot spin via bias voltage, gate voltage and an external applied 
magnetic field. The so prepared dot spin can be measured by it's imprint on 
the conductance of the quantum dot spin valve device as shown in Sec.0] We 
summarize our findings then in Sec. El 

2 Model Hamiltonian 

We describe the quantum-dot spin valve by the following Hamiltonian |12II13| : 

H = £ rk a cl ka C rka + Y £ n d l d n + U d\d ] d\d l (1) 

rka n 

+ H ( V rk an cl ka d n + /l.C. J . 
rkan 

The first term in Eq. JTJ treats the ferromagnetic leads (r = L/R) as large 
reservoirs of itinerant electrons. The Fermion creation and annihilation opera- 
tors of the lead r are labeled by , where k labels the momentum and a = ± 
the spin. The spin-quantization axis for the electrons in reservoir r is chosen 
along its magnetization direction m r . In the spirit of the Stoner model, the 
property of ferromagnetism is incorporated by assuming an asymmetry in the 
density of states £ a for majority (+) and minority (— ) spins. The degree of spin 
polarization in lead r is characterized by the ratio p r = (£r+ — £r— )/ (£r++£r-)- 
The lead magnetization directions ihl and mg can enclose an arbitrary an- 
gle 4>- Furthermore, the leads shall be so large, that the electrons can always 
be described as in equilibrium, i.e. with a Fermi distribution / r ((*>)■ An ap- 
plied bias voltage is taken into account by a symmetric shift of the chemical 
potential in the left and right lead by ±eV/2. 

The quantum dot can be modeled as an Anderson impurity, where S n and 
d n are the Fermion creation and annihilation operators of the dot electrons 
with the spin n =|,|- The spin quantization axis of the dot is, in general, 
chosen to be different from the quantization axes of both the left and the 
right lead. If an external magnetic field is applied, the spin quantization axis 
is chosen parallel to this field. The energy e„ of the atomic-like electronic level 
is measured relative to the equilibrium Fermi energy of the leads, and double 
occupation of the dot costs the charging energy U 3> k-gT. 

Electron tunneling between the leads and the dot is described by the last 
term in Eq. IjTJl. As we have chosen different spin quantization axes for the 
lead subsystems, parallel to the respective magnetization, the tunneling ma- 
trix elements V rk an are not diagonal in spin space. However, we require that 
tunneling is spin conserving. The tunneling amplitudes can, then, be sepa- 
rated in Vrkan = trk x U^ n , i.e. a spin-independent tunnel amplitude t r k and 
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a SU(2) rotation matrix V r an - The explicit shape of the matrix is determined 
by the choice of the dot spin quantization axis of the dot system. 

The tunnel-coupling strength is characterized by the transition rates 
r ra (uj) = 27r^ fc \t rk \ 2 S(uj — s r ka)- For simplicity, we assume the density of 
states £q, and the tunneling amplitudes t r to be independent of energy, which 
implies constant tunneling rates r ra . The spin asymmetry in the density of 
states in the leads yields spin-dependent tunneling rates, which are related to 
the leads' spin polarization by p r = (P r + — r r -)/(r r + + r r ~). 

Throughout this article we focus on the limit of dot-lead coupling r ra <C 
kBT,eV when transport is dominated by first-order tunneling, i.e., it is suf- 
ficient to calculate all expressions for the charge and spin current up to first 
order in r. This excludes the regimes of second-order transport (cotunneling) 
in the Coulomb-blockade region [2H| or the Kondo regime (see e.g. Ref. |%7|). 

3 Quantum-Dot-Spin Dynamics 

The tunnel coupling of the quantum-dot levels to spin-polarized leads yields 
a transfer of angular momentum across each of the tunnel junctions. This, 
together with a change of angular momentum due to an externally applied 
magnetic field, defines the dynamics of the quantum-dot spin polarization. The 
stationary value of the latter is determined by balancing all currents of angular 
momentum. As we discuss in detail below, the total spin current, i.e., transfer 
of angular momentum from the leads to the dot, consists of two qualitatively 
different contributions. One is associated with the fact that charge currents 
from or to a ferromagnet is spin polarized. This current, thus, transfers angular 
momentum along the magnetization directions of ferromagnet or quantum- 
dot spin accumulation. There is, however, also an additional contribution of 
transfer of perpendicular angular momentum, which can be expressed in terms 
of a many-body exchange field acting on the quantum-dot electrons. 

For a careful treatment of the total transfer of angular momentum, we first 
present a rigid calculation of the spin current through a single tunnel junction 
|28| in terms of non-equilibrium Keldysh Green's functions. This will help 
us to identify under which circumstances spin currents with perpendicular 
components will contribute. Afterwards, we specify our result to the weak- 
coupling regime of a quantum-dot spin valve and derive in this limit Bloch-likc 
rate equations for the quantum-dot spin. 

The currents will be functions of the unknown density matrix elements. 
To derive the stationary density matrix in a non-equilibrium situation, we set 
the change of charge and spin on the dot equal zero. Then, the Bloch equation 
for the spin and the continuity equation for the charge degree of freedom from 
a system of master equations. 
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3.1 Spin Current Expressed by means of Green's Functions 

Our calculation of the spin current in this subsection will be in close analogy 
to the derivation of the charge current according to Meir and Wingreen |29j . 
Let us first consider the spin current through one, say the left, tunnel barrier. 
For a clearer notation, we mostly drop the lead index in this section. The spin 
current Jl = (Jl) from the lead into the dot is defined by the negative of the 
time derivative of the total lead spin Sl = (ft/2) J2 ka p c\. a o ' a/3 c k/3 , where cr al3 
denotes the vector of Pauli matrices. From the Hciscnberg equation we get 



J L = --S L = -i 

dt ih 



(2) 



Making use of the Fermion commutation relation, we can find the spin-current 
operator as 

= 2J E Vkan.O'aP c kfj ~ Vkan eT af3 c fc/3 ■ (3) 
ka{3n 

By introducing the Keldysh Green's functions G^ k/3 (t) = i(c k p(0) dn(t)) , 
we can write the expectation value of the spin current as 

Jl = \ J2 J 2^ (y ka „(T* a pG< k p(u) - V kan cr Q p G^g n (o;)^ . (4) 

ka0n 

Since the Green's functions obey the Dyson equations 

f}< _ V V f n* C < —n < f 1 * 1 and = V* V* \ n < — 

^ka,n l—im ' ka,m\.yka^ J 'm,n i> ka u in ,n J all<a u n ,ka L-im ' ka,mVvka^' 'n,m 

9 ka Gn m ] i we can replace the Green's functions in Eq. with the dot Green's 
functions G„ m (t) = i(dj„ d n (t)) and the free Green's functions of the lead. The 
latter are given by g< a = 2irif+ (uj)S(ui - e ka ), g ka = -2irif^(u})5(uj - e ka ), 
gf* = l/( w - e ka + i0+), and gfj = Here, /+ stands for the Fermi 

distribution function in the lead L and =1 — f£. 

If we choose the dot spin quantization axis parallel to the lead magnetiza- 
tion we can substitute the tunnel matrix elements by Vk a , n = t k 5 an - After a 
lengthy but straightforward calculation, the spin current can be written as 



j l = i E / £ ^n n (r m + r n ) g> m + / L - ( W ) G 

m.n 



< 



"" r " VnnC^»n ^n) 



1 f lrn G n m (E) 



/i» (Gr m + G^) + - dE 



in E — to 



> (5) 



with the tunnel rates r n (v) = 2nJ2k l^rfcl 2 *^ — £rka)S an - 

This is the most general expression for the spin current flowing through a 
tunnel barrier. Since the Green's functions G n>m were not specified during the 
calculation, Eq. (J5J holds also for other electronic systems than single-level 
quantum dots. 
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If the dot state is rotationally symmetric about hil, all dot Green's 
functions G aa , non-diagonal in spin space vanish. Only in this case, the 
spin current is proportional to the difference between charge current IJ = 
i(e//i) / dui If [/l (w) G^ +f^(cu) G^. ] carried by spin- up electrons and charge 

current 1^ carried by spin-down electrons, 

J L = J£e, = ^(il-lt) (6) 

If the dot system breaks this rotational symmetry, for example due to spin 
accumulation along an axis different from hil, the simple result of Eq. © 
is no longer correct. In such a situation, the second line in Eq. JSJ yields an 
additional spin-current component, oriented transversal to both, the magne- 
tization of the lead, and the polarization of the dot. This spin-current compo- 
nent describes the exchange coupling between lead and dot spin, causing both 
to process around each other. Since the lead magnetization is pinned usually, 
only the dot spin precesses like in a magnetic field. 

Brataas et al. showed, that at normal metal - ferromagnet interfaces, 
incoming electrons, with a spin orientation non-collincar to the magnetization 
direction, may experience a rotation of the spin direction during backscatter- 
ing. This mechanism is described by the so called spin-mixing conductance, 
and also generates a transverse component of the spin current. 



3.2 Spin Current Between Ferromagnetic Lead and Quantum Dot 



We now specify the above expressions for a quantum-dot spin valve for weak 
tunnel coupling. Since the expression for the spin current in Eq. © does 
already explicitly depend linearly on the tunnel coupling P rai we only need 
the zeroth-order Keldysh Green's functions of the dot system to describe the 
weak-coupling regime. They are given by 



G^ 



-2mP s 5{u -e-U)- 2-kiP q 8{w - e) 
2-K\P a S{u) - e) + 2iriP d S(cu -e-U) 
2mP°6(u; -e - U) 
2iriP^S(tu -e) 

p<7 pa 



uj — e - 



i0+ 



U) — £ - 



(G a tny 



(7) 
(8) 
(9) 
(10) 

(11) 



where P* are the matrix elements of the reduced density matrix of the dot 
system, 



Pdot — 



/Fo \ 
' F T P} ' 



I 

P t L P l 



(12) 



V o o o pj 
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The diagonal, real entries P x = P* are the probabilities to find the dot in 
the state empty (0), occupied with a spin up (f) or down (J.) electron, or 
double occupied (d) with a spin singlet. The zeros in Eq. I|12[) in the off 
diagonals are a consequence of the total-particlc-number conservation. The 
inner 2x2 matrix is the SU(2) representation of the dot spin. The reduced 
density matrix contains five independent parameters. For convenience, we 
describe the quantum dot state by the probabilities for the three charge states 

T I 

Po, -Pi = Pj +Pp an d Pi (with the normalization condition Po+Pi+Pj = 1), 
and the average-spin vector S = (Pj + PpiP/ - iP/, P] - P/)/2. 

Similarly to deriving the expression for the spin current, we can get a 
formula for the charge current I r through tunnel contact r as 



I r = --^ / dw 

kan 



(VrkanG< rka (uj) -V* kan G< ka n (u)^ . (13) 



After choosing a spin-quantization axis for the dot spin and making use of 
the Dyson equation for the Green's functions, we can plug in the dot Green's 
function given above to obtain the result 



-Pr [f-(e) + f+(e + U)] S 



(14) 



that, of course, is independent of the choice of the dot spin's quantization 
axis. Here, we defined P r = (P r -f + P r jJ/2. 

It is worth to mention, that the dot spin S influences the conductance via 
the scalar product (S • m r ). Therefore the tunnel magnctorcsistancc depends 
cosine like on the relative angle enclosed by lead magnetization and spin po- 
larization, i.e. it just resembles the behavior of a tunnel junction between two 
ferromagnetic leads |%1 l?H 

The first-order spin current, on the other hand, is given by evaluating 
Eq. (JSJ), which leads to 

^ T S - pl(m r ■ S)m r 
j r = —l rPr m r — — — - + S x B r . (15) 

2e Tc,r 

The first term describes spin injection from the ferromagnetic lead into the 
quantum dot by a spin-polarized charge current. The injected spin is propor- 
tional to the lead polarization and the electrical current crossing the junction. 
This spin current contribution vanishes for vanishing bias voltage. 

The second term describes relaxation of the dot spin due to coupling to 
the leads. Since neither an empty nor a doubly-occupied dot can bear a net 
spin, the spin relaxation time r^ 1 = r r /h(l — / r (e) + f r (e + U)) equals the 
life time of the single-occupation dot state. This relaxation term is anisotropic 
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. The spin polarization of the lead suppresses the relaxation of a dot spin, 
which is aligned parallel to the lead magnetization. 

The third term in Eq. I|15|l describes transfer of angular momentum per- 
pendicular to the spin-polarization directions of lead and dot. The structure 
of this terms suggests the interpretation of B r as being an effective magnetic 
field that acts on the quantum-dot spin S. Its value, in the absence of an 
external magnetic field, is given by |12l \',V2\ 

B ,. = Pr ^f4^L 7 + £M). (M) 

■Kh J \uj — e — U uj - e J 

where the prime at the integral indicates Cauchy's principal value. From 
Eq. (|16J) it is clear that this field is an exchange field that arises due to the 
fact that the quantum dot levels are tunnel coupled to a spin-polarized lead. 
It is a many-body effect as all degrees of freedom in the leads contribute to the 
integral, and Coulomb interaction in the dot is important not to cancel the 
first with the second term in the integrand. The exchange field persists also 
for vanishing bias voltage. A signature of this exchange field in the Kondo- 
resonance splitting of transport through a single molecule has been observed 
recently |20] > with reported values of the field of up to 70 Tesla. 



3.3 Dynamics of the Quantum-Dot Spin 

We use the expressions for the charge and spin current, Eqs. Ijl4|l and (|15J) . to 
calculate the dynamics of the dot's charge and spin. Strictly speaking, above 
calculation holds only for static systems. To emphasize the physical origin of 
the following equations, we keep all time derivatives in this section, even if 
they should have the numerical value of zero. 

The continuity equation of the average dot charge (n) = ^ n nP n is given 

by 

Moreover, not only the total charge current through the dot is conserved, but 
also the charge current through the individual charge levels. Therefore we 
can split the charge continuity Eq. El hito the two contributions associated 
with transport processes in which either double occupied or an empty dot is 
involved. The affiliation to either contribution is indicated by the arguments of 
the Fermi functions, where the presence of the interaction energy U indicates 
processes with double occupation and the absence signals processes involving 
an empty dot. We get 



dPo 

dt 



^dT = E frie + U)Pi/2 - f-(e + U)P d - Pr f+(e + U)S ■ m r 



(18) 
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Similar to the charge continuity equation, the continuity equation for the 
dot spin reads 



dS 

~dt =Jl 



2e 



Jr. + S x B cxt 

Trel 

ES — p 2 (m r ■ S)m r 
l r p r m r 



S x B t 



S 

T"rel 



(19) 



with B to t = (Bl + Br + B cxt ). In addition to the spin currents entering the 
quantum dot from the left and right lead, there is one term describing spin 
precession due to an external magnetic field B cxt . It enters the equation in 
the same way as the exchange field with the left and right reservoir, so that all 
three of them add up to the total field B to t- Furthermore, we phcnomenologi- 
cally took into account the possibility of intrinsic spin relaxation, e.g., due to 
spin-orbit coupling, hypcrfine interaction with nuclei in the quantum dot, or 
higher-order tunnel processes such as spin-flip cotunneling, with a time scale 
T Te \. The total spin-decoherence time of the dot spin is, therefore, given by 



(t s ) 1 = (r rcl ) 1 + (t CjL ) 1 + (t c ,r)" 



(20) 




b) 

2 T 



Jl 











~~ s — 





— iH, — 

Fig. 2. a) The dot charge changes according to the electrical current through the 
tunnel barriers, b) The dot spin changes according to the spin currents through the 
tunnel barriers. In addition, an external magnetic field acts as additional source and 
the intrinsic spin relaxation as sink of angular momentum. 



The different handles to manipulate the quantum-dot spin are comprised 
the total field B to t- It contains the external magnetic field B ext as a direct tool 
to initiate a spin precession. However, also the the exchange fields can be used 
for this task ^3]- As we see from Eq. H16[). the exchange field depends on both 
the gate and bias voltage via the level position e and the Fermi distribution 
functions f r {u). 
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4 Manipulation and Detection of the Quantum-Dot Spin 
via Electrical Transport 

Since the spin state of the quantum dot enters the expressions for the charge 
current in Eq. 1)14(1. any manipulation on the quantum-dot spin can be de- 
tected in measuring the dc— charge current through the device. 

In order to calculate the charge current, we need to determine the sta- 
tionary solution for the density matrix, i.e. for the dot spin S and the charge 
occupation probabilities Pi. For these six variables, we need six independent 
equations: the probability normalization condition ^2 n P n = 1 ; the Bloch 
equation dS/dt = (contains three equations), and the two equations origi- 
nating from the charge continuity 

In the following we analyze stationary transport situations, i.e., neither the 
charge nor the spin of the dot changes with time, and the currents through 
the left and right tunnel junction are equal II = —Ir = I- For simplicity, we 
choose symmetric coupling Jl = Ir = -T/2, equal spin polarizations pl = 
Pr = p, and a symmetrically applied bias Vr = — Vl, = V/2 for the following 
discussion. 

In the following three subsection we consider the effect of the gate and 
transport voltage as well as an external magnetic field on the quantum-dot 
spin, and how this is reflected by electric transport through the quantum- 
dot spin valve. Then, in reversal, by experimentally measuring the transport 
characteristics of the device, one can conclude the spin state of the quantum 
dot. 

4.1 Gate Voltage Effect in Linear Response Regime 

To study the effect of the gate voltage on the quantum-dot spin via the 
gate- voltage dependence of the exchange field, we analyze the linear-response 
regime in the absence of an external magnetic field, and for simplicity with- 
out intrinsic spin relaxation. Without any applied bias voltage V = 0, i.e. 
in equilibrium the stationary solution of the rate equations for the charge 
occupation probabilities (|18|) is given by the Boltzmann distribution, P x ~ 
exp(— E x /k^T), and no current flows. Since the dot itself is non-magnetic, the 
dot spin vanishes S = 0. For a small bias voltage eV <C k^T, we can expand 
the master Eq. (|19fl and l|18|) up to linear order in V. With symmetric cou- 
plings to the left and right lead, the charge probabilities (Pq, Pi, Pa) become 
independent on V , thus, the occupation probabilities are given by their equi- 
librium value, but the spin degree of freedom is not. The linear charge current, 
which is polarized due to the lead magnetizations generate a finite dot-spin 
polarization along p(niL — hir), i.e., along the y axis in the coordinate system 
defined in Fig. |21 

The damping term in Eq. I|19|) limits the magnitude of spin accumulation. 
The term S x (B L + Br) yields an intrinsic precession of the spin around the 
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Fig. 3. Spin dynamics in the linear-response regime. Spin accumulates along the y 
direction. The spin precesses due to the exchange field that is along the x direction. 
Therefore, the stationary solution of the average spin on the dot is tilted away from 
the y axis by an angle a, plotted in Fig. 2[b). 



exchange field Bl +Br = Bo cos(</>/2) e x . In the steady state, the average dot 
spin is rotated by the angle 



Therefore the accumulated spin acquires both y and z components as seen in 
Fig. |21 This precession also reduces the magnitude of the accumulated spin to 



The precession angle a is plotted in Fig. EJb) as function of the level 
position e, that can be tuned by the gate voltage. The angle a changes its 
sign at e = —U/2, due to a sign change of the exchange field at this point. 
The level position e = — U/2 is special, since then the particle and hole like 
processes generating the exchange field compensate each other. 

As pointed out above, in the linear-response regime under consideration 
the charge occupation probabilities do not depend on the spin polarization 
of the leads. In particular, they are independent of the relative angle (ft of 
the leads' magnetization. This means that the eft dependence of the conduc- 
tance is determined by the product S • rriL = — S • itir, as can be seen from 
Eqs. (|14J) . It is the relative orientation of the accumulated spin and the drain 
(or source) that produces the (ft dependence of the current, rather than the 
product rriL • hir,, as in the case of a single magnetic tunnel junction. There- 
fore the (ft dependent linear conductance G lln = (dI/dV)\ v=o directly reflects 
the accumulated spin. The effect of the exchange field for the normalized 
conductance is seen from the analytic expression 




(21) 



| S | = pi t s cos a . 



(22) 



G lin ((ft) 
G lin (0) 




sin 2 (</>/2) 



(23) 



1 + OB t s ) 2 cos 2 W>/2) ' 
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Fig. 4. (a) Linear conductance normalized by r/ksT as a lunction of the level 
position e for different angles <j>. (b) Angle a enclosed by the accumulated spin and 
the y axis as defined in Fig. |3] (c) Derivation of the magnitude of the accumulated 
spin on the dot with respect to the source-drain voltage V . Further parameters are 
p = 0.9 and U = 10k B T. 



which is plotted in Fig. [5] for different values of the level position e. 

For e > 0, the quantum dot is predominantly empty, and for e + U < 
doubly occupied with a spin singlet. In this regions, the life time of a singly- 
occupied dot t c is short, and so is the lifetime of the dot spin. Therefore the 
rotation angle a is small and the normalized conductance as a function of the 
relative angle <f> of the lead magnetizations shows a harmonic behavior, see, 
e.g., the curve for e = 5fcsT in Fig.JSJ 

For —U < e < the dot is primarily singly occupied, so the spin dwell time 
is increased and the exchange field becomes important. It causes the above 
described spin precession, which decreases the product S-niL since the relative 
angle between rriL and S is increased and the magnitude of S is reduced. Thus, 
the spin precession makes the spin- valve effect less pronounced, leading to a 
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value of the conductance that exceeds the expectations made by Slonczcwski 
in Re. jH] for a single magnetic tunnel junction. 

For parallel and antiparallcl aligned lead magnetizations, = and <fi = ir, 
the accumulated spin and the exchange field also get aligned. In this case, the 
spin precession stops, even though the exchange field is still present. The <p- 
dependent conductance is not affected by the exchange field at this alignment, 
see Fig. |SJ 

4.2 Bias Voltage Effect in Non-Linear Regime 

We now turn to the non-linear response regime, eV > k^T, in order to discuss 
the effect of the bias- voltage dependence of the exchange field on the quantum- 
dot spin. Again, we assume that there is no external magnetic field, and no 
spin relaxation. 

In Fig. Efa) we show the current J as a function of the bias voltage V for 
an antiparallcl configuration of the leads' magnetizations and different values 
of the leads' spin polarization p. 

For non-magnetic leads, the current- voltage characteristic shows the usual 
Coulomb staircase. At low bias voltage, the dot is empty and transport is 
blocked. With increasing bias voltage, first single and then double occupancy 
of the dot is possible, which opens first one and then two transport channels. 
A finite spin polarization p leads to spin accumulation and, thus, to a reduc- 
tion of transport. A reduction of transport with increasing p is also seen for 
noncollinear magnetization. But there is a qualitative difference as can be seen 
in Fig. Etb) ■ A very pronounced negative differential conductance evolves out 
of the middle plateau as p is increased. To understand the negative differential 
conductance we first neglect the exchange field and then, in a second step, 
analyze how the exchange field modifies the picture. 
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Fig. 6. Current- voltage characteristics for antiparallel (a) and perpendicular aligned 
(b) lead magnetizations. Further parameters are Jl = -Th = r/2, Pl = Pa = P, 
e = 10fc B T, and U = 30fc B T. 



At the intermediate bias voltages, the dot can only be empty or singly 
occupied, double occupation is forbidden. Therefore all electrons entering the 
dot through the left barrier find an empty dot. In this regime the current 
I = (er/h)P explicitly depends only on the probability to find the dot empty. 
The transport through the dot must be suppressed by charge accumulation 
on the dot, i,c, Pi — ► 1 and Pq — ► 0. 

The origin of this charge accumulation becomes clear from the relation of 
charge and spin 



S=p 



r L p i-P 

— P m L — m R 

1 R I 



(24) 



If, in the steady state, the dot is primarily occupied by one electron, this elec- 
tron has a spin, which is antiparallel aligned to the drain lead. Due to this 
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antiparallcl alignment, the tunneling rate to the drain lead is maximally sup- 
pressed, while the tunnel coupling to the source lead is not as much affected. 
When the rate to the drain lead is weak, but strong to the source lead, then 
the dot is primarily occupied by one electron. 

The transport is suppressed, since an electron is trapped due to its spin 
alignment, and no second electron can enter the dot because of the Coulomb 
interaction. So this mechanism is a type of spin blockade but with a different 
physical origin compared to the systems described in literature jSHl EH ■ The 
suppression defines the local minimum of the current in Fig. EJb). At this 
point, the relevant exchange field component generated by the coupling to the 
left lead vanishes, so that spin precession becomes insignificant. Away from 
this point spin precession sets in as illustrated in Fig. The spin rotates 
about rriL, the spin blockade gets lifted, and the electron can now more easily 
leave the dot via the drain electrode. 



Fig. 7. For electrons polarized antiparallel to the drain lead, the influence of the 
effective field generated by the source lead is dominating. By rotating the spins, the 
spin blockade is lifted and therefore the conductance recovers. 



The particular value of the non-linear conductance is a consequence of 
the two competing effects. Spin blockade reduces, while spin precession again 
increases the conductance. Since the strength of the exchange field varies as 
a function of the level position with respect to the Fermi level, sec Fig.[Sfa), 
this recovery is non-monotonous, what leads to a negative differential con- 
ductance. To illustrate this further, we plot in Fig.|SJb) the current which we 
obtained when the spin precession contribution is in an artificial way dropped 
in Eq. (|19|) . and compare it with the total current. In the absence of the 
exchange field, a wide plateau is recovered, whose height is similar to the cur- 
rent, one would expect, if the lead magnetizations were aligned antiparallel. 
The peak at the left end of the plateau indicates that, once the dot level is 
close to the Fermi level of the source electrode, the spin blockade is relaxed 
since the dot electrons have the possibility to leave to the left side. 




xA 



B 
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Fig. 8. Panel (a) The absolute value of the effective exchange field contributions 
from the left and right leads. Panel (b) the current voltage dependence, with and 
without the influence of the exchange field. For both plots the parameters (j> = t/2, 
r L =r R = T/2, e = 10fc B T, U = 30fc B T, and p = 0.95 were chosen. 



However, this negative differential conductance occur only at relatively 
high values of the lead polarization. For symmetric tunnel coupling a spin 
polarization of p w 0.77 is needed, while for a strong asymmetry in the tunnel 
coupling the required spin polarization is reduced. 

The effect of the spin blockade on the ^-dependent of the current is de- 
picted in Fig- El We choose the bias voltage according to eV/2 = e + U/2, such 
that the influence of the exchange field is absent. For p = 0.5 still a sin 2 ^ 
dependence can be recognized. For higher values of the spin polarization the 
conductance drops faster and stays nearly constant at its minimal value due 
to spin blockade. This is just the opposite behavior than predicted for the 
linear-response regime as seen in Fig. [5] 

If such a high bias voltage is applied, that the dot can also be double 
occupied, the step-like behavior of the current voltage characteristic is recov- 
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Fig. 9. Angular dependence of the conductance with an applied voltage of V = 
e + U/2, i.e., the voltage generating the smallest influence of the exchange field. 
Further plot parameters are Fl = Fr, = F/2, ph = pr = p, e = 10/cbF, and 
U = 30fc B F. 



ered, see Fig. Elb). Away from the step, all appearing Fermi functions can 
be approximated by or 1, and following Eq. Ijl4(l the current is given by 
/ = (er/2h) [1 — pS ■ (niL — hir)]. Far away form the resonance, where the 
exchange field can be neglected, the accumulated spin is S — p(mj J — rriR)/4 
from which we get 

The suppression of transport due to the spin polarization p of the leads is 
comparable with the case of a single-tunnel junctions, when charging effects 
are of no importance. 

Wc close this section with the remark that while we plotted only results for 
the case e > 0, in the opposite case e < the current-voltage characteristics 
is qualitatively the same. 

4.3 External Magnetic Field 

In the previous subsections were studied quantum-dot spin dynamics that 
is evoked by the exchange field. But one also make use of an externally- 
applied magnetic field B cxt . It turns out that with an external field one can 
measure the spin-decoherence time T2. To emphasize this point, we explicitly 
allow intrinsic spin relaxation on the dot. Then, we observe a separation of 
the charge life time r" 1 = r"^ + and the spin life time on the dot 
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An external field leads to the Hanlc effect JJ, i.e., the decrease of spin 
accumulation in the quantum dot due to precession about a static magnetic 
field. Indeed, this was the effect used by Johnson and Silsbee [35] and others 
|17j to prove non-equilibrium spin accumulation. 

Optical realizations of such Hanlc experiments |36| always involve an en- 
semble averaging over different dot realizations, so the outcome of the mea- 
surement is TJ rather than Ti- By measuring the Hanle signal via the conduc- 
tance through a quantum dot attached to ferromagnetic leads, this ensemble 
averaging is avoided. 

In a recent experiment Zhang et. al. |25j realized this kind of setup but 
with a whole layer of aluminum dots in a tunnel junction between two Co 
electrodes. Even so the measurements involve averaging over different realiza- 
tions of the dots, multi levels and local magnetizations, they clearly observe 
a Hanle resonance in the magnetoresistance of the device. 

For simplicity we assume symmetric couplings Pl = Pr, equal degree of 
lead polarizations j»l = PR = P and consider the linear-response regime only. 
There is a variety of possible relative orientations of the external field and the 
leads' magnetizations to each other. In the following, we consider two specific 
cases in detail, as they are convenient to extract useful information about the 
spin-decoherence time in one case, and to prove the existence of the exchange 
field in the other one. 



Antiparallel Aligned Lead' Magnetizations 

We first focus on two ferromagnetic leads with magnetization directions anti- 
parallel to each other, see Fig. an d an arbitrary aligned external field. The 
configuration has the advantage that the exchange field contributions from 
the two leads cancel, and the spin dynamics is only govern by the external 
field B cxt . The linear conductance, then, is 

Go T C 1 + (B oxt T s )2 ' 1 > 

where Go = e 2 Pi/r c k^T is the asymptotic value of the conductance for a large 
magnetic field, |B oxt | — > oo, for which the spin accumulation is completely 
destroyed. The latter is proportional to the single occupation probability Pi. 

If we assume the field to be aligned perpendicular to the lead magne- 
tizations (see Fig. I10[) . we find the Lorentzian dependence on the external 
magnetic field that is familiar from the optical Hanle effect. The depth of the 
dip is given by p 2 t s /t c while the width of the dip in Fig. 1101 provides a direct 
access to the spin lifetime r s . Of course, the conversion of applied magnetic 
field to frequency requires the knowledge of the Landc factor g. which must 
be determined separately like in Ref. |22| . 
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Fig. 10. Differential conductance, for ferromagnetic leads with anti-parallel mag- 
netization, as a function of the magnetic field to applied perpendicular to the accu- 
mulated spin. The half line width of the Hanle resonance directly determines the 
spin-decoherence time t s . 



Magnetic Field Applied Along hil + niR 

Finally, we discuss the case of a non-collinear configuration of the leads' mag- 
netizations with a magnetic field applied along the direction hil + mjt as 
shown in Fig. 1111 

In this case, both the exchange field and the external magnetic field arc 
pointing along the same direction niL + mR, so their magnitude is just added. 
The linear conductance is, then, 



where <j> is the angle enclosed by niL and hir. The conductance as function of 
applied magnetic field as plotted in Fig. lllf a) reaches its minimal value when 
the sum of external and exchange field vanishes. The exchange field leads to a 
shift of the minimum's position relative to |B cx t| = |13| . In real experiments, 
depending on the particular sample geometry, one can expect a magnetic stray 
field, which is not considered to be part of the experimentally applied magnetic 
field B oxt . These stray fields also lead to a shift of the conductance minimum. 
However, the analyzed setup of external field and magnetizations' directions 
allows for a stringent experimental verification of spin precession due to the 
exchange field. To separate the exchange field from the influence of possible 
stray fields its gate voltage dependence can be used. The exchange interaction 
as function of the dot gate voltage is plotted in the inset of Fig. lTlT b'). While 
the stray fields does not depend on gate voltage, the exchange field does. In 




T C 1 + (B«t + B L + B R ) 2 T S 2 ' 



(27) 
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Fig. 11. Linear conductance of the dot for an applied external magnetic field B ox t 
along niL + hir. a) Linear conductance as a function of the applied field for e = 0. 
b) Linear conductance as a function of the level position e without external field 
(dotted) and for the applied field |B cx t| = O.lF/h (solid). Further parameters are 
4> — 37r/4, p — 0.8, U = 7k^T, and r re i = 0. The vertical lines relate the conductance 
increase of the dot at e = for a magnetic field Hujb = O.ir. 



the flat band limit it even changes sign as a function of gate voltage. By 
plotting the conductance as function of the gate voltage in Fig. Illf b). we 
can observe the typical Coulomb blockade oscillations, when the energy level 
of the empty or singly-occupied dot becomes resonant with the lead Fermi 
energy The interplay of exchange and external field leads to an increase of 
conductance for one resonance peak, but to a decrease for the other resonance. 

5 Conclusions 

We discussed the possibility to generate, manipulate, and probe single spins 
in single-level quantum dots coupled to ferromagnetic leads. A finite spin- 
polarization of the quantum-dot electron is achieved by spin-polarized charge 
currents from or to the leads at finite bias voltage. Any manipulation of the ac- 
cumulated spin, e.g. by an external magnetic field or by an intrinsic exchange 
field, is detectable in the electric current through the device. The occurrence 
of the exchange field is a consequence of many-body correlations that are one 
of the intriguing features of nanostructures with large Coulomb interaction. 

We determine the dynamics of the quantum-dot spin by deriving expres- 
sions for the spin currents through the tunnel barriers. In addition to a con- 
tribution that is associated with the spin-polarization of the charge currents 
from or to ferromagnets, there is a second contribution describing transfer of 
angular momentum perpendicular to the leads' and dot's magnetization that 
can be expressed in terms of the exchange field. 
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In order to manipulate the quantum-dot spin we suggest to make use of the 
gate- and bias voltage dependence of the exchange field or to apply an external 
magnetic field. In particular, the spin precession modifies the dependence of 
the linear conductance on the opening angle of the lead magnetizations. The 
degree of modification is tunable by the gate voltage. In nonlinear response, 
the bias-voltage dependence of the exchange field can give rise to a negative 
differential conductance. An application of a tunable external magnetic field 
allows one to determine the dot-spin lifetime and to verify the existence of 
the intrinsic spin precession caused by the exchange coupling. 
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